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مقايسة عملکرد مدل‌های سری زمانی 54181314 و Holt-Winters‏ با روش‌های ھوش 
مصنوعی در بيش بينى طوفان‌های گردوغبار (مطالعة موردی: استان سیستان و بلوچستان) 
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جکیده 

هدف از این پژوهش مقايسة عملکرد مدل‌های سری زمانی Holt-Winters  SARIMA‏ با روش‌های هوشس 
مصنوعی شامل شبكة عصبی مبتنی بر توابع ab‏ شعاعی (RBF)‏ و سیستم استنباط عصبی - فازی تطبيقى 
(ANFIS)‏ به‌منظور پیش‌بینی فراوانی روزهای همراه با طوفان گردوغبار (FDSD)‏ در فصل آتی است. بدین منظورء 
از داده‌های ساعتی گردوغبار و کدهای سازمان جهانی هواشناسی در ينج ایستگاه سینوبتیک استان سیستان و 
بلوجستان با طول 5555 آماری ۲۵ساله (۲۰۱۴-۱۹۹۰) استفاده شد. نتایج نشان داد روش ۸۸۳15 نسبت به ساير 
روش‌هاء بهترین عملکرد را داشت و معیارهای ارزیابی MAE RMSE R‏ و NS‏ آن به‌ترتیب از ۰/۷۲ ۰/۵۷ 
۲۳ ۱ تا ۰/۹۵ ۰۰۱۵۱ 9۰۰/۴۰ ۰/۹۶ متغیر بود. همچنین. با افزایش شاخص متوسط FDSD‏ در 
ایستگاه‌ها (از ۱/۰۶ تا ۷/۱۱ دقت بيش بينى همة روش‌ها افزایش داست. بر همین اساس» در سری زمانی 
۸ ضريب همبستگی بين مقادیر مشاهداتی و پیش‌بینی‌شده شاخص 7959 از ۰/۶۴ به ۰/۷۹ افزایش 
یافت. برای روش‌های RBF dolt-Winters‏ و ANFIS‏ مقدار نیز ضریب همبستگی ias‏ تیب از ۰/۷۰ تا e [AV‏ 


۵۹ تا ۰/۹۲ و ۰/۷۲ تا ۰/۹۵ متغیر بود. درمجموع. با مقايسة روش‌های مورد استفاده. روش هوش مصنوعی 
5 بهترین و مدل‌های سری زمانی SARIMA‏ و Holt-Winters‏ بدترین عملکرد را داستند. 


وازگا نکلیدی: تحلیل سری زمانى» شبکه‌های عصبی, گردوغبار» ۸۷و PACF‏ نواحی بحران یگردوغباری. 


مه 


مشد مه 

يديدة گردوغبار یکی از مخاطره‌های طبیعی است و سالانه خسارت‌های زیادی در جنوب غرب و جنوب شرق کشور ایجاد 
می کند. گردوغبار یکی از منابع آلوده‌کنندة هوای پاک و عاری از آلودگی است و علت افزايش طوفان‌های گردوغبار در 
ایران به دلیل موقعيت جغرافیایی کشور با توجه به قرارگیری در کمربند خشک و نیمه‌خشک جهان است و از سوی ديكر 
نزدیک‌بودن با کشورهایی همچون عربستان» عراق» و سوریه که دارای سطح وسیعی بیابان به دلیل عامل‌های مختلفی 
است. مانند کمبود پوشش گیاهی. بارش «Sul‏ و دمای بالا (عبدالشاه‌نژاد و همکاران» ۱۳۹۹: (VY‏ همچنین» وقوع 
خشک‌سالی در یک منطقه با تأمين آب در أن منطقه dll)‏ تنگاتنگی دارد و هر جه اقليم منطقه خشک‌تر باشد» احتمال 


' نويسندة مسئول» تلفن: ۰۹۱۲۲۱۵۹۷۳۸ Email: AliaghatQut.ac.ir‏ 


به طور كلىء منابع عمدة گردوغبار عبارت‌اند از: سرزمينهاى كمارتفاع و کمبارانی كه از ميانكين بارش ساليانه كمتر 
از ۲۵۰ ميلىمتر برخوردارند. كشور ایران نيز از این قاعده مستثنى نيست (گودی» ۲۰۱۳: ۱۰۱). ايران» به ee‏ قرارگرفتن 
در كمربند خشك و نیمه‌خشک» در معرض سيستمهاى متعدد كردوغبار محلى و فرامنطقهاى است. همچنین» نيمة غربى 
کشور به دليل عرصه‌های بيابانى و بادهاى منطقه‌ای» هرساله أسيبهاى زيادى از اين مخاطره را تجربه م ىكند (حسين 
حمزه و همكاران: ۱۳۹۵: [AY‏ 

در بين متغیرهای اقلیمی» بارش و دما از اهمیت ب‌سزایی برخوردارند و شبیه‌سازی دقیق بارش و دما و ارتباط أن با 
گردوغبار در مطالعات هواشناسی و کشاورزی از اهمیت زیادی برخوردار است. با وجود شبیه‌سازی بارش و دما به دليل 
وجود همبستگی درونی بين آن‌ها مشکل است (تانارته و همکاران» ۲۰۱۲). نتایج علمی مشابه توسط هرویجر و همکاران 
(۲۰۰۷) و هاننبرگر و نیکول (۲۰۱۴) نشان داد با افزايش دما تعداد وقوع طوفان‌های گردوغبار نیز افزايش مىيابد. 
همچنین, نتایج اين تحقیق در ارتباط با همبستگی پارامترهای اقلیمی با طوفان‌های گردوغبار همسوست با نتایج مطالعات 
sj‏ و همکاران (۲۰۰۶) و آمگالان و همکاران (۲۰۱۷)» که افزایش دماء افزايش سرعت ob‏ و کاهش بارندگی را از 
مهم‌ترین دلایل افزايش رخداد پديدة گردوغبار معرفی کرده‌اند. آیگنگ و همکاران (۲۰۰۹) با بررسی ارتباط گرمایش 
جهانی و تغییرات طوفان‌های گردوغباری جينء به اين نتيجه رسیدند که گرمایش جهان سبب کاهش شيب عرضی دما و 
در نتيجه کاهش شدت باد می‌شود. 

جمالی‌زاده تاج‌آبادی و همکاران (۱۳۸۶) توانایی دو روش شبکه‌های عصبی مصنوعی (ANN)‏ و ماشین‌های بردار 
پشتیبان (SVM)‏ در پیش‌بینی طوفان‌های گردوغبار شهر زابل را با طول دورة آماری ۲۶ساله (۲۰۰۵-۱۹۸۰) بررسی 
کردند. از exl‏ كه استفاده از نوع 7-5۷ و تابع کرنل åh‏ شعاعی (RBF)‏ بهترین جواب‌ها را داده است» توصیه شده در 
مطالعات دیگر نیز قابلیت‌های اين روش‌ها بررسی شود. 

سبحانی و همکاران (۱۳۹۴) به بررسی گردوغبار و ارزیابی امکان پیش‌بینی أن بر اساس روش‌های آماری و مدل 
۳ در زابل با آمار چهارده‌ساله پرداختند. نتایج پیش‌بینی گردوغبار با مدل ANFIS‏ نشان‌دهندة cub‏ بالای أن در 
پیش‌بینی گردوغبار در اين ایستگاه است. علی‌یاری و همکاران (۲۰۰۸) نیز به پیش‌بینی کوتاه‌مدت آلودگی هوا به کمک 
شبکه‌های عصبی پرسپترون چندلایه J(MLP)‏ خط حافظه‌دار تأخير (DLM)‏ گاما و سیستم استنتاج فازی- عصبی 
تطبیقی (ANFIS)‏ پرداختند. نتيجه نشان داد روش پیشنهادی ترکیبی بر اساس بهینه‌سازی ازدحام ذرات (PSO)‏ و فیلتر 
کالمن برای آموزش ANFIS AS,‏ توانایی مناسبی در بهبود عملکرد پیش‌بینی دارد. 

سروستان و فلاح قالهرى (YA)‏ به پیش‌بینی يديدة گردوغبار در هشت ایستگاه منتخب استان خوزستان با استفاده 
از مدل‌های سری زمانی باکس جنکینز پرداختند. نتيجه نشان داد الگوهای مناسب ماهانه به‌ترتیب برای رامهرمز 
ARIMA (Y e ANAN)‏ آغاجاری (۰۱۰۱()۱۰۱۰۱ ARMA (Y‏ بهبهان ARIMA (ENFAN)‏ آبادان (Y V)‏ 
۸ دزفول ARIMA (MY) y)‏ امیدیه (۸۳(۸۰۸) «ARIMA‏ اهواز (۰۸۳(۸۸) ARIMA‏ و 
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Adaptive Neuro-Fuzzy Inference System 
. Multi Layer Perceptron 

. Delay Line Memory 

. Particle Swarm Optimization 


مقايسة عملکرد مدلهاى سرى زمانى SARIMA‏ و Holt-Winters‏ با روش‌های هوس مصنوعى در پیش‌بینی e‏ ۵۶۹ 


مسجدسليمان (۰۳۰۴()۱۰۱۰۱) ARIMA‏ می‌باشند كه اين مدل‌ها از دقت بسيار خوبى برای پیش‌بینی كردوغبار برخوردار 
نبودند و همجنين بيش بينى تعداد روزهاى طوفان گردوغبار برای سال‌های ۲۰۱۸ تا ۲۰۲۷ نشان داد كه «ele ll‏ آبادان 
و مسجدسليمان بيشتر با يديدة كردوغبار مواجه‌اند. 

گودرزی و روزبهانى (۱۳۹۶) كارايى مدل‌های سرى زمانی ARIMA‏ و هالت وينترز در پیش‌بینی دما و بارش ماهانة 
ایستگاه لتیان را بررسی کردند. gelo‏ پیش‌بینی توسط مدل‌های ARIMA‏ فصلی نشان از قابلیت بالای اين مدل‌ها در 
پیش‌بینی دمای ماهانه داشت و مدل ۲( ۰۰۱۲ (۰.۰۰۱) ARIMA‏ به‌عنوان مدل سری زمانی مناسب برای داده‌های دما 
تشخیص داده شد. اگرچه مدل هالت- وینترز نيز از دقت (e Yo‏ در پیش‌بینی دما برخوردار بوده نسبت به مدل ARIMA‏ 
خطای بیشتری داشت. همجنین» نتایج نشان داد مدل‌های ARIMA‏ فصلی و مدل هالت-وینترز قابلیت بالایی در 
پیش‌بینی مقادیر بالای بارش ماهانه در ایستگاه لتیان ندارند. اين مدل‌ها عمدتاً مقادیر بارش‌های پایین و متوسط را بهتر 
برآورد می کنند و تخمینگر مناسبی در بارش‌های بالا نیستند. 

سبحانی و همکاران (Y)‏ به بررسی مدل‌سازی و پیش‌بینی گردوغبار در غرب ايران پرداختند. بدین منظور از داده‌های 
گردوغبا do»‏ و رطوبت YA‏ ایستگاه مناطق درگیر شدید با گردوغبار در ايران در 85b‏ زمانی ٩۲ساله‏ (۲۰۱۸-۱۹۹۰) استفاده شد. 
بر اساس نتایج پژوهش. در بهترین شرایط مقدار RMSE‏ مدل ANFIS‏ برابر با ۱۱,۶۷ و مدل RBF‏ برابر با ۲,۱۹ است. بنابراین» 
قدرت RBF cà»‏ در پیش‌بینی گردوغبار در سال‌های شبیه‌سازی‌شده بیشتر است. بر اساس نتایج خروجی مدل RBF‏ در 
پیش‌بینی گردوغبار برای سال‌های آتی ایستگاه‌های مورد مطالعه. در هر دو مقیاس میانگین و حداکثر فراوانی كردوغبار, 
ایستگاه‌های غربی و جنوب غربی منطقة مورد مطالعه بیشتر در معرض گردوغبار در سال‌های أينده قرار گرفتند. 

انتظاری و سروستان (۱۳۹۶) به بررسی گردوغبار و پیش‌بینی أن در استان خوزستان برای باز زمانی ۲۰۲۳-۲۰۱۹ با 
a an IST‏ ی اساسا دای ار باه شام مات ارس 
خوزستان در طول دورة آماری ۲۰۱۰-۱۹۹۰ استفاده شد. نتیجة یافته‌های مناسب‌ترین مدل‌های سری زمانی گردوغبار 
برای شهرهای مسجدسلیمان و بهبهان مدل پیش‌بینی هالت-وينترزء شهر دزفول به‌تنهایی از مدل ARIMA (Y ۰۰ Y)‏ و 
شهرهای اهواز رامهرمز آغاجاری, آبادان» و اميديه دارای مدل نمو هموار ساده به‌دست آمد. 

تغییر اقلیم و بررسی اثرهای أن بر منابع آب از موضوعات مهمی به‌شمار می‌رود که پژوهشگران زیادی هدف تحقیق 
خود را در اين زمينه تعريف کرده‌اند. بنابراین» آهمیت پیش‌بینی وقوع پدیده‌های جوی نظیر كردوغبارء که بسیار حاثر 
همیت است. بیش از پیش حس می‌شود. تأثیر يديدة گردوغبار در ایران أن قدر وسیع است که بیش از نیمی از 
استان‌های کشور را به نحوی با مسائل و محدودیت‌های اين يديدة طبیعی درگیر کرده است که» علاوه بر اثرهای 
زیست‌محیطی» موجب اختلال در اجرای طرح‌های توسعة پایدار ملی شده و تاکنون پیامدهای منفی زیادی به دنبال 
داشته و خواهد داشت. از طرفی» بیشتر مطالعات داخلی در اين زمينه مربوط به روند پدیده‌های گردوغبار در محدوده‌های 
کوچک. مطالمات سینوپتیکی, و ماهوارای آن است. culla‏ با توجه به اينکه این پدیده اثرهای تامطلوب و پیامدهای 
منفی در زمینه‌های اجتماعی» اقتصادی» و سلامت مردم داشته است. مطالعه و پیش‌بینی و ارتباطسنحی أن با نوسانات 
اقلیمی ضروری است. در اين مطالعه به مقايسة عملکرد مدل‌های سری زمانی SARIMA‏ و Holt-Winters‏ با روش‌های 


هوش مصنوعی RBF‏ و ۸۱۷۳5 به‌منظور پیش‌بینی فراوانی روزهای همراه با طوفان گردوغبار (FDSD)‏ پرداخته شده 


1. Frequency of Dust Stromy Days 


است و بنابراين اين تحقيق می‌تواند راهنما و الكوى جديدى به‌منظور پیش‌بینی يديدة طوفان كردوغبار در كشور باشد؛ 
زيرا تا كنون چنین مطالعه‌ای با اهداف فوق‌الذ کر انجام نشده است. 


مواد و روس‌ها 
منطقة مورد مطالعه و داده‌های مورد استفاده 


در اين پژوهش به بررسی مقايسة عملکرد روش‌های آماری کلاسیک SARIMA)‏ و (Holt-Winters‏ و مدل‌های هوش 
مصنوعی RBF)‏ و (ANFIS‏ به‌منظور پیش‌بینی فراوانی روزهای همراه با طوفان گردوغبار (FDSD)‏ در ينج ایستگاه 
سینوپتیک استان سیستان و بلوچستان (زابل زاهدان, ايرانشهرء خاش, و سراوان» که بیشترین فراوانی روزهای همراه با 
طوفان‌های گردوغبار در کشور را شامل می‌شوند (عراقی‌نژاد و همکاران؛ ۱۳۹۷: ۲۱ با جامعة آماری ۲۵ساله (۱۹۹۰- 
۴ در مقیاس فصلی abu‏ شد. بدین منظور از داده‌های ساعتی قدرت دید افقی و کدهای سازمان جهانی 
هواشناسی (WMO)‏ استفاده شد. مشاهدات پدیده‌های هواشناسی به فاصلة سه ساعت یک بار و درمجموع هشت بار در 
شبانه‌روز ثبت می‌شود. در اين مشاهدات» پدیده‌های بصری أب و هوا بر اساس دستورالعمل سازمان جهانی هواشناسی در 
اک که TEENS E‏ کی SO E‏ رقن ET‏ کار کر انمق قفا 
مختلف هواشناسی از ۱۱ كد استفاده می‌شود (اولینگ‌سای و همکاران» ۲۰۱۴). کدهای مربوط به طوفان‌های گردوغبار در 
جدول اارائه شده است. با توجه به اينكه داده‌های گردوغبار به‌صورت تفکیک‌نشده از سایر پدیده‌های هواشناسی از 
سازمان هواشناسی در اختیار قرار می‌گیرد» در وهلة اول کدهای مربوط به گردوغبار از ساير پدیده‌های اقلیمی در ساعات 
مختلف همدیدی جدا شد (گودی و میدلتون. ۲۰۰۶). بر اساس تعریف سازمان جهانی هواشناسی روزهای همراه با 
طوفان‌های گردوغبار به روزی گفته می‌شود که حداقل در یکی از هشت سینوپ (گزارش‌های سه‌ساعتة ديدهبانى) یکی از 
کدهای مربوط به گردوغبار (۰۶» ۰۷ ۰۸ ۰۹ ۲۰ تا YO‏ و (AA‏ در بخش هوای حاضر گزارش شده باشد. البته, به شرط 
اينكه داده‌های قدرت دید افقی متناظر با آن كد گردوغبار کمتر از ۱۰۰۰ متر به ثبت رسیده باشد (گودی و میدلتون» 
EIE nua‏ شیرتا ماک رفن ایا و قرو ا شا T‏ 
هواشناسی گردوغبار استفاده شده است. 


جدول ۱. کدهای سازمان جهانی هواشناسی مرتبط با فرسایش بادی و پدیده‌های گردوغبار (اولینگ‌سای و همکاران. ۲۰۱۴) 


مه ناشی از گردوغبار 

id‏ غبار يا شن برخاسته از زمين 

۸ طوفان كردوغبار 

v‏ طوفان گردوغبار اتفاق‌افتاده در گذشته (اتفاق‌افتاده در يك ساعت قبل از مشاهده يا در ايستكاه) 

۳۰ گردوغبار خفيف يا متوسط شن و ماسه همراه b‏ كاهش ميدان ديد كمتر از ٠٠٠١‏ متر اما بيشتر از ۲۰۰ متر 

Y‏ گردوغبار پایدار يا متوسط شن و ماسه همراه با کاهش میدان دید کمتر از ٠٠٠١‏ متر اما بیشتر از ۲۰۰ متر 

۳۲ شروع يا افزايش گردوغبار خفيف يا متوسط شن و ماسه همراه b‏ كاهش ميدان ديد كمتر از ٠٠٠١‏ متر اما بیشتر از ٠٠١‏ متر 
۳۴ طوفان گردوغبار پایدار شدید همراه با کاهش دید کمتر از ۲۰۰ متر 

۵ شروع يا افزايش طوفان گردوغبار شدید همراه با کاهش دید کمتر از ۲۰۰ متر 


1. World Meteorological Organization 


مقايسة عملكرد مدلهاى سرى زمانى Holt-Winters  SARIMA‏ با روش‌های هوس مصنوعی در پیش‌بینی e‏ ۵۷۱ 


يس از انتخاب ایستگاه‌ها و بررسى داده‌ها در بازة زمانى ۲۵ساله (۲۰۱۴-۱۹۹۰)» تعداد روزهاى همراه با طوفان 
كردوغبار (FDSD)‏ براى ينج ایستگاه هواشناسى مورد مطالعه در استان سيستان و بلوچستان با استفاده از داده‌های ديد 
افقى و كدهاى سازمان هواشناسی محاسبه و در جدول Y‏ نشان داده شده است. ضمن اينكه ایستگاه‌های هواشناسی. 
طول و عرض جغرافيايى» ارتفاع از سطح درياء متوسط شاخص FDSD‏ در مقیاس فصلی و تعداد روزهای گردوغباری 
به‌ترتیب نزولی قابل مشاهده است. ایستگاه زابل» که در بين ابستگاه‌های مورد مطالعه در مرزی‌ترین نقطة استان سیستان 
و بلوچستان واقع شده با ۷۱۱ روز در بازة زمانی ۲۵ساله به‌عنوان رکورددار بیشترین روزهای همراه با طوفان گردوغبار در 
اين استان ثبت شده است. بر اين اساس» ایستگاه‌های خاش 9 سراوان به‌ثر تیب b‏ ۱/۸۰ 9 ۱۰۶ روز کمترین فراوانی 
روزهاى همراه با طوفان گردوغبار را در بين ابستگاه‌های مورد مطالعه به خود اختصاص دادهاند. بر اساس اين توضیحات 
به صورت جغرافیایی می‌توان بیان کرد که با پیشروی از مناطق شرقی و شمالی استان (مناطق مرزی) به سمت جنوب و 
جنوب غرب استان از وقوع فراوانی روزهای همراه با طوفان گردوغبار کاسته شده است. ede‏ أن را می‌توان در وزش 
بادهای ۱۲۰روزه» كه در نتيجة تشکیل سیستم پرفشار در شمال شرق دریای خزر و شکل‌گیری یک سامانة کم‌فشار در 
خشكشدن تالابهاى منطقه جستوجو كرد. نتايج این بخش b‏ مطالعات مهرابى و همكاران (WAY)‏ انتظاری و 
سروستان (YAT)‏ عراقی‌نژاد و همکاران (۱۳۹۷). و انصاری قوجقار و همکاران (۱۳۹۹) همخوانی دارد. شكل ۱ موقعيت 


جدول ۲. فراوانی روزهای همراه با طوفان‌های گردوغبار استان سیستان و بلوچستان (۲۰۱۴-۱۹۹۰) 
عرض جغرافیایی طول جغرافیایی ميانكين فصلى FDSD FDSD‏ 


ce‏ (درجه) ibd au ME‏ روز روز 
y ۱۹۵ 2-7 yyy- EUR‏ ۱.۶ 
خاش Me vis WAY BAM YAY‏ 
Yvy PEOR‏ ۲ع ۵۹۱ YM YA‏ 
زاهدان ۳۹۲۳۸ ۶۰۵۳ .۱۳۷ YA‏ ۳۴ 
زابل YM ۴۳۸۹ ۹ yyy-‏ ۷/۱ 
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آناليز اوليه داده‌ها 

برای أناليز» مدلسازىء و پیش‌بینی از نرمافزارهاى R‏ و Matlab‏ استفاده شد. پیش از هرگونه اقدام در جهت مدل سازى يا 
پیش‌بینی سرى زمانی» تجزيه و تحليل توصيفى اولیه‌ای از داده‌ها امرى ضرورى است. در اين مطالعه أناليز اولية داده‌ها بر 
روی مجموعه داده‌های اموزشی صورت گرفت که شامل شناسایی ویژگی‌های مهمی نظیر همبستگی.» الگوهای فصلی. 
تغییرات دوره‌ای تودرتو در الگوهای فصلىء روند» داده‌های پرت» و هرگونه نوسانات قابل توجه دیگر در اين سری‌هاست. 
هه ddl‏ ادوا bs‏ ررد کو کور از ا مستا نون انش ات که cnl usb‏ اسان اما چ مات فا گید 3 
قاری ترش di a‏ باق io odo E So odia o‏ 
پیش‌بینی فرض می‌کنند داده‌های سرى زمانى ايستا می‌باشند؛ زيرا با این فرض يبجيدكى رياضى مدل‌های برازش 
داده‌شده كاهش می‌یابد؛ حال أنكه اغلب در واقعيت چنین نيست. در اين مطالعه از آزمون رياضى Dickey-Fuller‏ برای 
تشخیص ایستایی gum b‏ دادهها استفاده شد (YY : ۰۰۵۵ «u$ sS)‏ آنگاه در مواردى كه CS y»‏ زمانى ناايستا بود از 
تفاضل گیری " و تبدیلات توان برای حذف ناایستایی روند يا الگوهای فصلی استفاده شد. آنالیز اولي داده‌ها نیز از طریق 
تجزیه و تحلیل بصری نمودارهای زمانی شامل ACF‏ ۳۸۲۳ و محاسبات آمار توصیفی پایه انحام گرفت. برای 
شناسایی الگوهای موجود در سری‌های زمانی از تابع Decompose‏ در نرم‌افزار R‏ استفاده شد که اين تابع از روش 
میانگین متحرک برای تحزیة سری زمانی به اجزای فصلی» روند» و بخش نامنظم و تصادفی استفاده می کند. 


انتخاب و برازس مدل‌ها 

دارد. تحقیقات گسترده‌ای انحام شده است که نشان می‌دهند یک روش پیش‌بینی برتر برای همه موارد و در حالت goo‏ 
پیش‌بینی با صحت قابل قبول» بايد روش‌های پیش‌بینی مختلفی مقایسه شود و توانایی هر یک در پیش‌بینی مشاهدات 
نمونه‌های آینده ارزیابی شود. به دلیل وجود الگوهای فصلی در داده‌های ایستگاه‌های مورد مطالعه» مدل‌های پر کاربرد و 
دارای صحت بالاتر در مطالعات مشابه. كه قابلیت بيش بينى الگوهای فصلی را نیز داشتند. به‌عنوان مدل‌های پیش‌بینی 
سرى زمانى در اين مطالعه درنظر گرفته شده‌اند كه اين مدل‌های آماری عبارت‌اند از مدل‌های سرى زمانى Holt-‏ 
Winters‏ و SARIMA?‏ (ابوكىسارفو و همکاران» ۲۰۱۵) و مدلهاى يادكيرى ماشين شامل دو روش RBF‏ و ANFIS‏ 


(سينك و ميشراء ۲۰۱۵). 


مدل SARIMA‏ 
مدل SARIMA‏ گسترش يافتة مدل ARIMA?‏ برای داده‌های فصلى است. منظور از سری‌های زمانى فصلی» سرىهايى 
اهوت كه 2525 ous‏ تم Delos jio‏ ات obs‏ وا قرام o‏ که کی کی کا قمعا من سرض 
زمانی» كه داراى الكوهاى فصلئاندء نهتنها بررسى وجود همبستگی بين مشاهدات فعلى و مشاهدات اخير ضرورى است 


]. Stationary 

2. Differencing 

3. Auto Correlation Function 

4. Partial Autocorrelation Function 

5. Seasonal Autoregressive Integrated Moving Average 
6. Average Autoregressive Integrated Moving 
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که بایستی وجود همبستكى ميان مشاهدات فعلى و مشاهدات قبلى در فصلهاى گذشته نيز بررسى شود كه روشى براى 
شناسايى وجود يا فقدان الگوی فصلى در داده‌ها نيز هست. اين مدل فرض می‌کند که سرى زمانى درنظ osiad S‏ خطى 
است و از یک توزیم آماری شناخته‌شدة خاص, مانند توزیع نرمال» تبعیت می‌کند. محبوبیت مدل SARIMA‏ عمدتاً به 
دلیل انعطاف‌پذیری أن در نمايش ساده انواع مختلفی از سری‌های زمانی برای فرایند ساخت مدل بهینه است. برای اجرای 
مدل SARIMA‏ از الگوریتم ارائه‌شده توسط هایندمن و خندکار استفاده شده است (هایندمن و خندکان ۲۰۰۸). 


Holt-Winters مدل‎ 

۰ * e 6 ۰ 5 ^ 3 n 
اين مدل یکی از محبوبترين روش‌های هموارسازى نمايى است. اين مدل بر اساس الكوهاى فصلى و روند موجود در‎ 
داده‌ها ایحاد مىشود. محبوبيت اين روش به دليل مزايايى از قبيل نياز محاسباتی كمترء سادكى استفاده» و دقت پیش‌بینی‎ 
(Multiplicative و‎ Additive) این مدل به‌خصوص برای سرىهاى دارای الگوی فصلى است. اين روش دو كونه معادله‎ 
استفاده شد که برای سری‎ Additive دارد. از آنجاکه نوسانات سری‌های زمانی ابن مطالعه تسيا ثابت بردتت از معاد‎ 
انس‎ B) به صورت‎ lx: t= l,..., n! زمانى‎ 
X "ua Zl, + hb, + E (Y) 


كه h‏ افق پیش‌بینی است. برای مثال» در اين مطالعه كه پیش‌بینی برای فصل آتی aw)‏ ماه آینده) است» h‏ برابر با Y‏ بوده 
و ساير پارامترها عبارت‌اند از: 

Level: /, za(y, — s, ,,) ۲ (d-a) I, , b, 4) 

Trend: b, =B* (I, — La) +A -Bb (Y) 

Seasonality : s, 2yCy, — 1, 4— b, 4j) - (1—y)s, 4 


كه B* ax‏ و y‏ پارامترهای هموارسازی‌اند و ja‏ به صورت 1 > © > 0 1 > *8 > 0, و © -1 > ۲ > 0 بوده و 
همجنين m‏ طول 8555 فصلى است. 


شبكة عصبى مبتنی بر توابع پایة شعاعی (RBF)‏ 

شبكة عصبى مبتنى بر توابع ياية شعاعی (RBF)‏ داراى پاية رياضياتى بسيار قوی بر مبنای فرضية منظم‌سازی برای حل 
مسائل مشكل است. به طور كلى اين شبكه از سه لايه شامل à Y‏ ورودی» مخفی» و خروجى تشكيل شده است (چن و 
همکاران. ۱۹۱: ۲۰۳). در اين شبکه تابع انتقال كوسى در لابة مخفى استفاده می‌شود و در ALY‏ خروجى تابع انتقال 
خطى است. RBF‏ يايه و اساس نوع خاصى از شبکه‌های عصبى به نام شبکه‌های عصبى آماری است. در شكل Y‏ یک 
شبکة RBF‏ نشان داده شده است. همان طوركه در شكل مشحص iul‏ نورون RBF‏ یک تابع گوسی FT‏ ورودى اين 
تابح فاصلة آقلیدسی بين هر ورودی به نورون با بردار مشخ صشدة هم‌اندازة بردار ورودی است. اين تابح گوسی از رابطة 


f(X,b)-e' r-|X, - X,|*0.8326/h (V) 


]. Hyndman 
2. Khandakar 
3. Exponential Smoothing 


در اين رابطه X,‏ ورودی شبكه با خروجى نامعلوم» ,ل ورودی مشاهداتى در زمان يا مكان sib‏ / پارامتری است كه 
پهنای تابع كوسى را كنترل می کند. خروجى اين تابع بين صفر تا یک متغير است. محاسب۵ خروجى ,۲ بر اساس متغير 


مستقل X,‏ به صورت زیر به‌دست مى أيد: 
Y. = LW * f(X „,b) + Bias (Y)‏ 
در این رابطه LW‏ و Bias‏ بهترتيب وزن ماتريس ارتباطات بين AY‏ مخفى و LY‏ خروجى و ماتريس باياس AY‏ 


rac c op 


E f CX, $ D) = e Vx — x, || 0.8326 AF 
| 227 | 
X3, ل‎ RE 
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سيستم استنباط عصبى - فازى تطبيقى (ANFIS)‏ 
نمع استباط aa‏ فازی تطبیقی ba (ANEIS)‏ نكن ور سال ۱۹۹۴ معرفی شد ANEIS 1 scil)‏ 


شبیه به یک شبكة عصبی چندلابه است؛ با اين تفاوت که علاوه بر الگوریتم‌های بادگیری شبکة عصبی. از منطق فازی 
نیز بهره می‌گیرد. یک مدل ANFIS‏ از ينج لايه تشکیل شده است؛ اين ينج لابه به‌ترتیب عبارت‌اند از: AY‏ ورود 
اطلاعات à Y‏ محاسبة وزن قوانین فازی. LY‏ نرمال‌سازی وزن‌های قوانین به‌دست‌آمده, Y‏ محاسبة قوانین, LY‏ 
جمع‌بندی و خروجی شبکه. در اين تحقیق تابح عضویت درنظرگرفته‌شده تابح عضویت دوزنقه‌ای و الگوریتم آموزش 
شبکه روش هیبریدی درنظر گرفته شد. شکل ۲ شماتیکی از ANFIS‏ را نشان می‌دهد. 





شکل Y‏ شماتیک سیستم استنباط عصبی- فازی تطبیقی (ANFIS)‏ 
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مدل‌های پیش‌بینی 

به‌منظور پیش‌بینی کمی فراوانی روزهای همراه با طوفان‌های گردوغبار (FDSD)‏ چهار مدل مختلف استفاده شد. سپس» 
بر اساس معیارهای ارزیابی بهترین مدل در هر بيش بينى انتخاب شد. متغیر تعداد فصل (های) گذشته برای پیش‌بینی 
فصل آتی به‌کار گرفته شد که می‌تواند شامل یک دوه سهء و چهار فصل قیل از متفیر فراوانی روزهای همراه ب 
طوفان‌های گردوغبارباشد. در aa‏ برای پیش‌بینی متفیر فراوانی روزهای همراه با طوفان‌های گردوغبار در فصل بعده 
از سری زمانی ۳5۲ با یک گام تأخیر استفاده شد. بدین صورت که برای پیش‌بینی در زمان 1+ از مقدار آن در زمان t‏ 
استفاده شد. در مدل ۲ پیش‌بینی در فصل بعد بر اساس مقادیر فراوانی روزهای همراه با طوفان‌های گردوغبار تا دو فصل 
قبل استفاده شد و به طور مشابه در مدل ۳ و ۴ پیش‌بینی‌ها بر اساس مقادیر متفیر مورد نظر تا سه و چهار ab‏ قبل 
صورت پذیرفت. روابط ۸-۵ زیر در تکمیل توضیحات بیان‌شده ارائه می‌شود. شکل ۶ مدل‌های مورداستفاده در انتخاب 
فصل (های) گذشته برای پیش‌بینی در فصل آتی را نشان می‌دهد. روش‌های بيش بينى به کاررفته در اين پژوهش شامل 
RBF Holt-Winters ۸‏ و ANFIS‏ است. 





FDSD, ب‎ f (FDSD,,) ١لدم‎ (à) 
FDSD, „=f (FDSD,,,, FDSD +) مدل۲‎ )۶( 
FDSD, „=f (FDSD,,, FDSD, ,,, FDSD 4) ۲ مدل‎ (Y) 
FDSD, „=f (FDSD,,, FDSD,_, FDSD, ,رو‎ FDSD,, 4) مدل۴‎ (A) 
"o 
(Y) مدل‎ 
UNIS SARIMA 
Holt-Winters 
(f) مدل‎ RBF 
© ANFIS 

سكل ۴. ساختارهاى مورد استفاده در پیش‌بینی كردوغبار 

آماده‌سازی داده‌ها 


يس از بررسی و کنترل کیفی آمار ایستگاه‌ها و رفع نواقص آماری» همگنی داده‌ها با استفاده از آزمون Run Test‏ ارزیابی 
و تصادفی‌بودن داده‌ها در سطح اطمینان ۹۵درصد پذیرفته شد. در ادامه» با توجه به وجود سری زمانی متغیر فراوانی 
روزهای همراه با طوفان‌های گردوغبار اقدام به پیش‌بینی أن با استفاده از روش‌های .Holt-Winters .SARIMA‏ 
RBF‏ و ANFIS‏ به‌منظور استفاده از روش‌های هوش مصنوعی شد. از ۸۰ درصد داده‌ها برای مرحلة آموزش و ۲۰ 


در صد برای آزمون استفاده شد. پارامترهای مناسب رو شهاى مختلف از طریق ازمون و خطا تعیین شد. همچنین. 


بدمنظور جلوكيرى از بیش‌برازش از صحت‌سنجی متقاطع استفاده شد. در روش صحت‌سنجی متقاطع در هر بار اجرای 
مدل برای پیش‌بینی با یک پارامتر مشخصء یک يا تعدادی از داده‌های مجموعه آموزش كنار گذاشته می‌شوند (به‌عنوان 
داده‌های صحت‌سنحی) و مدل بر اساس داده‌های باقی‌مانده آموزش می‌بیند. سپس داده‌های کنارگذاشته‌شده به مجموعة 
آموزش برگردانده شده و داده يا مجموعة دیگری از داده‌های آموزش به‌عنوان دادة صحت‌سنحی از اين مسير جدا می‌شود 
و بار دیگر شبکه آموزش داده می‌شود. اين فرایند تا زمانی که همه داده‌های آموزش در قسمت صحت‌سنحی استفاده 
شوند» ادامه می‌یابد. يس از أن بر اساس میانگین خطای مدل در تخمین داده‌های صحت‌سنحی, پارامتر (های) مناسب 
شبکه عصبی تعیین می‌شوند. برای پیش‌بینی از چهار مدل معرفی‌شده فوق الذکر استفاده شد. مدل مناسب بر اساس 
معیارهای ارزیابی برای هر روش در مقیاس فصلی استخراج شد. پس از «ol‏ تحلیل نتایج بر پایه بهترین مدل برای هر 
يك از روش‌ها صورت پذیرفت. 


معیارهای ارزیابی 


e. 1 ١ eus ۰ ۰ ۰ ۰ * مه $ همه ^ ا‎ m 
جذر ميانكين‎ (R) برای ارزيابى مدلها و مقايسة نتايج روش‌های مختلف از جهار معيار ارزیابی يعنى ضريب همبستكى‎ 
استفاده شد. اين معيارها به‌صورت‎ (NS) و نش‌ساتکلیف‎ (MAE) ميانكين قدر مطلق خطا"‎ (RMSE) مربعات خطا"‎ 


زیر تعريف می‌سوند: 





در اين رابطه‌هاء 01 polio‏ مشاهداتی در گام زمانی eli‏ 11 مقادیر پیش‌بینی‌شده در گام زمانی لام 0 ميانكين مقادیر 
مشاهداتى» Í‏ ميانكين مقادیر پیش‌بینی‌شده» و n‏ تعداد داده‌هاست. هر مدل که دارای R‏ و NS‏ بیشتر و RMSE‏ و 


شكل ۵ clus‏ شاخص FDSD‏ را در پنج ايستكاه سينويتيك منتحب استان سیستان 9 بلوچستان از منظر فراوانى 
روزهای همراه با طوفان گردوغبار نشان می‌دهد. بر اساس اين «Is‏ شاخص FDSD‏ در ایستگاه‌های سراوان» خاش» 


1. Correlation 

2. Root Mean Square Error 
3. Mean Absolute Error 

4. Nash-Sutcliffe 


مقايسة عملكرد مدلهاى سرى زمانى SARIMA‏ و Holt-Winters‏ با روش‌های هوس مصنوعی در پیش‌بینی AYY e‏ 


افزايش تعداد روزهای گردوغباری در ایستگاه زابل از پراکندگی نوسانات کاسته دة 9 بك شدت ]4 ن افزوده شده انش 9 
مقادير اوج گردوغبار به صورت متمرکز در كنار یکدیگر مشاهده می‌شوند که حکایت از وقوع طوفان‌های گردوغبار 
بی‌دریی در اين ایستگاه از سال ۲۰۰۰ به بعد دارد. با توجه به شکل ۵ وجود فصول متوالی عاری از طوفان‌های 
گردوغبار انقطاع و داشتن یک سری زمانی منفصل, و به تبع أن تشکیل طوفان‌هایی با فاصلة زمانی بیشتر نسبت به 
ایستگاه زابل در سال‌های ۳ تا ۲۰۰۱ در ایستگاه‌های سراوان» خاش ایرانشه 9 زاهدان مشهود sa‏ بر اساس 
۲۰۱۲-۲۰۸ را در ایستگاه‌های مورد مطالعه تشخیص داد. نتایج مشايهى از ویژگی‌های شاخص FDSD‏ در مقیاس زمانی 
فصلی توسط زینالی (۱۳۹۵) عراقی‌نژاد و انصاری قوجقار (۱۳۹۶) و عراقی‌نژاد و همکاران (۱۳۹۷) ارائه شده است. 
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كام زمانی (فصلی) 


سکل ۵. سری زمانی ساخص FDSD‏ در مقیاس زمانی فصلی برای ایستگاه‌های مورد مطالعه 


اولين قدم در پیش‌بینی سری‌های زمانی» شناسایی ویژگی‌های مهمی نظیر همبستگی, الگوهای فصلی. تغییرات 
دوره‌ای در الگوهای فصلی. روند و ابستایی يا غیرایستایی داده‌هاست. وجود همبستگی در ميان داده‌های سری زمانی» 
امکان اجرای مدل‌سازی و پیش‌بینی مقادیر بعدی سری زمانی را فراهم می‌کند. در اين مطالعه. وجود همبستگی با رسم 
نمودارهای ACF‏ و PACF‏ برای هر ایستگاه با فاصلة اطمینان ۹۵ درصد (خطوط نقطه‌چین) بررسی شد (شکل ۶). اگر 
معنی‌داری با صفر اختلاف دارد. همان‌طور که در نمودارهای ACF‏ و PACF‏ ایستگاه‌های مورد مطالعه مشاهده می‌شود. 
وقفه‌های زمانی معنی‌دار : نشان‌دهندة وجود همبستگی بين مقادیر زمانی است كه اجرای مدل‌سازی 9 پیش ‌بینی مقادیر 
أينده (فصل آتی) شاخص ۴۵52 را برای هر ينج ایستگاه مورد مطالعه ممکن می‌سازد. گام مهم در تجزیه و تحلیل 
سری‌های زمانی بررسی ایستایی يا غیرایستایی میانگین سری‌های زمانی است. با توجه به شکل ۶ نمودارهای ACF‏ و 
PACF‏ ایستگاه‌ها حکابت از ایستایی داده‌های همه سری‌های زمانی به غير از ابستگاه زابل دارد. 
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Partial Autocorrelation Function for saravan 
(with 5% significance limits for the partial autocorrelations) 








Partial Autocorrelation Function for khash 
(with 5% significance limits for the partial autocorrelations) 
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Partial Autocorrelation Function for iranshahr 
(with 5% significance limits for the partial autocorrelations) 
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Partial Autocorrelation Function for zahedan 
(with 5% significance limits for the partial autocorrelations) 
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Partial Autocorrelation Function for zabol 
(with 5% significance limits for the partial autocorrelations) 
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Partial Autocorrelation 
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پژوهش‌های جغرافياى طبیعی» 5,93 ۵۲ شمارة F‏ زمستان ۱۳۹۹ 


Autocorrelation Function for saravan 
(with 596 significance limits for the autocorrelations) 
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(with 5% significance limits for the autocorrelations) 
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Autocorrelation Function for iranshahr 
(with 5% significance limits for the autocorrelations) 
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Autocorrelation Function for ۵۱ 
(with 5% significance limits for the autocorrelations) 
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Autocorrelation Function for zabol 
(with 596 significance limits for the autocorrelations) 
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شكل ۶. نمودار خودهمبستگی (ACF)‏ و خودهمبستگی جزئی (PACF)‏ شاخص FDSD‏ در ایستگاه‌های مورد مطالعه 


















































مقايسة عملكرد مدلهاى سرى زمانى Holt-Winters  SARIMA‏ با روش‌های هوس مصنوعی در پیش‌بینی e‏ ۵۹ 


در داده‌های غيرفصلى و سالانه 00-75 با يك بار تفاضل‌گیری درجة أول داده‌ها ايستا می‌شوند؛ ولی» با توجه به 
فصلی‌بودن و مقیاس داده‌هاء با تفاضل‌گیری درجة دوم (D22)‏ داده‌های ایستگاه زابل ایستایی شدند که نمودار توابع 
خودهمبستگی (ACF)‏ و خودهمبستگی جزتی (PACF)‏ أن بعد از تفاضل‌گیری در شکل V‏ نشان داده شده است. با توجه 
به توابع خودهمبستگی و خودهمبستگی جزئی» محدودة تغییرات آتورگرسیون و ميانكين متحرک تعیین و با استفاده از 
معیارهای ارزیابی مناسب بهترین مدل سری زمانی برای هر ایستگاه استخراج شد. همچنین, از آزمون Dickey-Fuller‏ 
نيز برای تشخیص ایستایی يا ناایستایی سری‌های زمانی بر روی مجموعه داده‌های آموزشی استفاده شد (جدول Y‏ در 
آزمون 1016-۳1167 سطح معنی‌داری برای احتمال ۰,۰۵ P-Value‏ درنظر گرفته شد. فرض این آزمون نالیستایی 
سری زمانی را بررسی می کند. بر اساس اين آزمون» صرفاً سری زمانی ایستگاه‌های زابل ناایستاست که نتايج نمودارهای 
ACF‏ و PACF‏ ایستگاه‌های مورد مطالعه را تأييد می کند. 


Autocorrelation Function for zabol.new Partial Autocorrelation Function for zabol.new 
(with 5% significance limits for the autocorrelations) (with 596 significance limits for the partial autocorrelations) 


Autocorrelation 
Partial Autocorrelation 











شكل V‏ نمودار خودهميستكى (ACF)‏ 9 خودهمیستکی جزئنی (PACF)‏ بعد از تفاض ل كيرى فصلى مرتيه دوم (D-2)‏ در ایستگاه زابل 


جدول .Y‏ نتايج آزمون Dickey-Fuller‏ 


ابستكاه P-Value‏ 
سراوان JA‏ 
خاش ا 
ایرانشهر ۱ 
زاهدان 7“ 
زابل yA‏ 


جدول ۴ polio‏ معیارهای ارزیابی مربوط به مرحلة آزمایش روش‌های مختلف به‌منظور پیش‌بینی فراوانی روزهای 
همراه با طوفان گردوغبار را نشان می‌دهد. نتایج پیش‌بینی شاخص 170510 حاکی از عملکرد مناسب روش‌های سیستم 
استنباط عصبی- فازی تطبیقی (ANFIS)‏ و شبكة عصبی مبتنی بر توابع ياية شعاعی (RBF)‏ در همة ایستگاه‌های مورد 
مطالعه (زابل» «ual;‏ ایرانشهر» خاش و سراوان) است. بر اساس جدول ۴ در ایستگاه‌های زابل» زاهدان» و ایرانشهر 
روش‌های سیستم استنباط عصبی- فازی تطبیقی (ANFIS)‏ شبكة عصبی مبتنی بر توابع Alb‏ شعاعی (RBF)‏ مدل‌های 
سری زمانی dHolt-Winters‏ و SARIMA‏ به‌ترتیب رتبۂ اول تا چهارم بهترین روش‌های بيش بينى شاخص FDSD‏ را بر 
اساس معیارهای ارزیابی به خود اختصاص دادند. برخلاف ایستگاه‌های مذكورء در ایستگاه خاش بعد از روش ۸۱۳۲5 


وو 














استفاده از مدل سرى زمانى ۲101-۷۷10675 برای بيش بينى شاخص FDSD‏ عملكرد بهتری نسبت به روش ۸8۴ از خود 
نشان داد. به طورى كه ضرايب همبستگی (R)‏ و نش ساتكليف (NS)‏ بين polio‏ مشاهداتى فراوانى روزهاى همراه با 
طوفان كردوغبار و مقادير پیش‌بینی‌شدة أن بدترتيب از ۰۷۱ به ۰۷۴ و ۰۷۰ به ۰,۷۱ و همجنين معیارهای خطاى 
RMSE‏ و MAE‏ بدترتيب از ۰۸۷۰ به ۰,۶۹ و ۰,۵۶ به ۰,۵۵ بهبود يافت. همچنین» در ايستكاه سراوان نيز بعد از روش 
5 روش RBF‏ و مدل سرى زمانى Holt-Winters‏ عملكرد تقريباً يكسانى از خود نشان دادند؛ به طورى که 
معيارهاى ارزیابی MAE RMSE R‏ و NS‏ در روش RBF‏ نسبت به روش Holt-Winters‏ بدترتيب از ۰,۶۹ ego‏ 


وو 


«RÀ 9 ۵‏ به uS -7Y 9 ۵۶ A e.‏ يافت. b‏ نوجه بك اينكه معيارهاى ارزيابى RMSE‏ و NS‏ نشان دهندة 
برتری روش RBF‏ نسبت به مدل Holt-Winters‏ و معیارهای ارزیابی R۸‏ و MAE‏ حكايت از برتری مدل Holt-Winters‏ 
ایستگاه داشت. علت بهبود نتایج مدل سری زمانی Holt-Winters‏ در ایستگاه‌های خاش و سراوان را می‌توان در کاهش 
وقوع فراوانی روزهای همراه با طوفان گردوغبار جست‌وجو کرد؛ به طوری که اين موضوع حاکی از قابلیت بالای مدل 
سری زمانی Holt-Winters‏ در بيش بينى مقادیر کم شاخص FDSD‏ است؛ (Ja‏ با توجه به آهمیت موضوع پیش‌بینی 
طوفان‌های گردوغبار در مناطق بحرانی کشور (پرتعداد» نمی‌توان چندان به توانایی‌های اين مدل در پیش‌بینی فراوانی 
روزهای همراه با طوفان گردوغبار اتکا کرد. همجنينء نتایج نشان داد که مدل سری زمانی SARIMA‏ نسبت به ساير 
روش‌ها قابلیت (e Yo‏ در پیش‌بینی شاخص FDSD‏ در هیچ یک از ایستگاه‌های مورد مطالعه ندارد. نتایج اين بخش با 


جدول ۴. معیارهای ارزیابی محاسبه‌شدة مدل‌های آماری کلاسیک سری زمانی و روش‌های هوش مصنوعی برای پیش‌بینی شاخص FDSD‏ 


حالت معیارهای ارزیابی سراوان خاش ایرانشهر ‏ زاهدان زابل 
مدل بهينه Y F ۳ ۴ F‏ 
NA (Y) ۶۱ (Y) ۰,۷۰ (Y) ۶۹ (Y) /۶۶ (Y) R‏ 
(Y) «NA )۴( «FA (Y) NY (Y) ۷۵ (Y) RMSE SARIMA‏ ۰۶۵ 
(Y) ۰۵۷ (Y) ۰,۵۶ (Y) ۰,۵۷ (Y) ۰,۶۰ (Y) MAE‏ ۰۸۵۵ 
(Y) ۰,۶۰ (Y) ۰,۶۸ (Y) ۰,۶۵ (Y) ۰,۶۲ (Y) NS‏ ۰۸۱ 
مدل بهینه ۳ ۴ ۳ ۳ F‏ 
JAY (Y) ۰,۶۲ (Y) ۸۷۵ (Y) NS (Y) ۸۷۰ (V) R‏ 
(Y) ۰۸۷۸ (V) ۰,۶۷ )۲( ۰,۶۹ (Y) ۰۸۷۲ )۲( RMSE Holt-Winters‏ ۰,۶۲ 
(Y) ۰۸۵۲ (Y) ۰۵۲ (Y) ۰۸۵۵ (Y) ۰,۵۶ (Y) MAE‏ ۰۸۵۱ 
(V) ۰,۶۲ (v) ۰۸۷۰ (v) ۰۸۷۱ (Y) ۰,۶۷ (N) NS‏ ۰۸۵ 
مدل بهينه Y y ۳ ۳ Y‏ 
(Y) ۰۶۷ (Y) ۰,۷۶ (Y) ۰,۷۱ (Y) ۰۶۹ (Y) R‏ ۰۹۲ 
(Y) ۰۸۷۳ (Y) ۰۶۵ (Y) ۰,۷۰ (Y) ۰,۷۱ (Y) RMSE RBF‏ ۰۶۰ 
(Y) ۰,۵۰ (Y) ۰,۴۹ (Y) ۰,۵۶ (Y) ۰/۵۹ (Y) MAE‏ ۰,۴۸ 
(Y) ۰,۶۵ (Y) ۰۸۷۱ (Y) ۰۸۷۰ )۳( ۰,۶۸ (Y) NS‏ ۰/۱ 
مدل بهينه Y Y Y Y Y‏ 
(Y) ۰۸۷۸ (Y) ۰۸۷۷ (5) ۰۸۷۶ (V) ۸۷۲ (V) R‏ ۰/۵ 
(Y) ۰۸۵۲ )۱( ۰۵۲ (Y) ۰,۵۳ (Y) ۵۷ )۱( RMSE ANFIS‏ ۰,۸۵۱ 
(Y) ۰,۴۱ (Y) ۸۴۱ (Y) ۰,۴۱ (Y) FY (Y) MAE‏ ۰,۴۰ 
UNE )۱( ۰,۷۵ (V) NY (Y) ۷۲ (Y) ۰,۷۱ )۱( NS‏ 
متوسط روزهاى همراه با طوفان گردوغبار در مقياس فصلى ۱,۰۶ ۸ vM vM YM‏ 
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با توجه به اينكه در جدول Y‏ ترتیب ایستگاه‌ها به لحاظ شاخص FDSD‏ از ایستگاه سراوان تا زابل به صورت صعودی 
odg‏ و متوسط شاخص FDSD‏ در مقیاس فصلی از ۱,۰۶ تا ۷,۸۱۱ متغیر است» می‌توان دریافت که عملکرد هم روش‌های 
پیش‌بینی شاخص ۳5۲ با افزايش تعداد روزهای همراه با طوفان گردوغبار رابطة مستقیمی دارد. بدین معنی که با 
افزايش فراوانی روزهای همراه با طوفان گردوغباره مقادیر همه معیارهای ارزیابی مورد استفاده بهبود می‌یابند؛ به نحوی 
که در روش SARIMA‏ ضریب همبستگی بين polio‏ مشاهداتی و بيش بينىشدة شاخص FDSD‏ از ۰/۶۶ به ۰۷۹ 
افزايش می‌یابد. همچنین» برای روش‌های RBF Holt-Winters‏ و ANFIS‏ نيز به‌ترتیب از ۰۷۰ تا CUN‏ ۰,۶۹ تا ۰/۲ 
و ۰/۷۲ تا ۰/۹۵ در ایستگاه‌های مورد مطالعه متغیر است. نتيجة دیگری که از جدول Y‏ قابل استنتاج است عملکرد بهتر و 
مناسب همة روش‌های پیش‌بینی شاخص ۳۲5۲ بر اساس معیارهای ارزیابی در ایستگاه ایرانشهر نسبت به ایستگاه 
زاهدان است؛ اين در حالی است که فراوانی روزهای همراه با طوفان‌های گردوغبار در ایستگاه ایرانشهر نسبت به ایستگاه 
زاهدان کمتر است. علت اين موضوع را می‌توان در سری کامل و بدون انقطاع شاخص FDSD‏ در ایستگاه ایرانشهر 
نسبت به ایستگاه زاهدان جست‌وجو کرد؛ به طوری که در ایستگاه ایرانشهر تنها در دو مقطع زمانی تابستان و پاییز 
سال‌های ۲۰۰۵ و ۲۰۱۱ دو فصل متوالی بدون روزهای همراه با طوفان گردوغبار رخ داده است. 

در جدول Y‏ اعداد داخل پرانتز روبه‌روی polio‏ معیارهای ارزیابی رتبه‌بندی روش‌های مختلف پیش‌بینی فراوانی روزهای 
همراه با طوفان‌های گردوغبار را نشان می‌دهد. بر این due‏ در همة ایستگاه‌های مورد مطالعه به‌منظور پیش‌بینی شاخص 
۴02 روش سیستم استنباط عصبى-فازى تطبیقی (ANFIS)‏ بیشترین تعداد رتبة اول (بهترین حالت پیش‌بینی) را داشته 
اسان یاههار :رز راتت E E EE E CUTS‏ يكل 
وجود خواهد داشت. بدين ترتیب مدل‌های سرى زمانى dHolt-Winters SARIMA‏ و روش ۸8۴ اصلاً رتبة نخستی به 
خود اختصاص نداده‌اند؛ در حالى كه روش ANFIS‏ ۲۰ بار در رتبة نخست قرار گرفت. درمحموع و با اين حساب» روش 
5 با اختلاف زیادی بهتر از روش‌های دیگر پیش‌بینی بوده و بالاترین دقت پیش‌بینی را داشته است. نکتة دیگری که 
بايد به أن اشاره کرد عملکرد مناسب شبكة عصبی مبتنی بر توابع Alb‏ شعاعی (RBF)‏ است. هرچند که اين روش در دو 
ایستگاه خاش و سراوان به‌ترتیب عملکرد ضعیف‌تر و lg‏ یکسانی با مدل سری زمانی Holt-Winters‏ داشته در غیاب 
روش 4۸۴15 با W‏ بار قرارگیری در رتبة دوم» بهترین عملکرد را به خود اختصاص داد. همچنین, نتایج دیگری از جدول ۴ 
قابل برداشت است که انتخاب مدل‌های ييجيدهتر به‌عنوان مدل بهينة پیش‌بینی کننده در ایستگاه‌های مورد مطالعه است؛ به 
نحوی که برای بيش بينى شاخص FDSD‏ در هم ایستگاه‌ها از مدل ۳ و L) Y‏ سه و چهار گام تأخیر) استفاده شده است. با 
توجه به اینکه بادهای ۱۲۰روزه عامل اصلی وقوع طوفان‌های گردوغبار در شرق و جنوب شرق کشورند» افزایش سرعت و 
E‏ یهام مزر رن tona‏ و ad dioe sus‏ درا ویک هرد 
قدرتمند برای وقوع طوفان‌های بعدی بوده كه در مطالعات فرج‌زاده اصل و علیزاده (۱۳۹۰)» راشکی و همکاران 4Y* W)‏ 
سبحانی و همکاران (Y)‏ یارمرادی و همکاران (۱۳۹۷)» عراقی‌نژاد و همکاران (WAY)‏ و قربانی و مدرس (WA)‏ بدان 
اشاره شده است. اين ذرات به‌جامانده از طوفان‌های قبلی می‌تواند دليل مهمی برای تأثير طوفان‌های گردوغباری چند فصل 
گذشته برای شکل‌گیری طوفان‌های گردوغبار در فصل آتی باشد. 

شکل ۸ polio‏ مشاهداتی و پیش‌بینی شاخص FDSD‏ را برای داده‌های آزمایش نشان می‌دهد. اين شکل بپبود 
عملکرد مدل‌های مورد استفاده به‌منظور پیش‌بینی فراوانی روزهای همراه با طوفان گردوغبار را با به کارگیری روش 
5 در همه ایستگاه‌ها به‌خوبی بیان می کند. به طوری که در همة ایستگاه‌های سینویتیک مورد مطالعه در استان 


سيستان و بلوجستان» روش ANFIS‏ بهترین مطابقت بين مقادير پیش‌بینی و مشاهداتى شاخص FDSD‏ را ull»‏ 
به ۴۵ درجه (مدلهاى سرى زمانى SARIMA‏ و Holt-Winters‏ و روش (RBF‏ به یک خط دقيقاً ۵ درجه (روش 
(پیشروی از ایستگاه سراوان do‏ سمت ايستكاه زابل)» ارتباط بين مقادیر مشاهداتی 9 پیش‌بینی‌شده در همه روشها 


(مدل‌های سری زمانی و روش‌های هوش مصنوعی) مطابقت بیشتری با نيمساز ربع اول بيدا می‌کند. 


| dua 
US 


۴1 مقدار 


كه 


n 
T 


ز تبس ببس 
فا 
Cc‏ 


^ 





4 10 1 
1 j 4 L1 ANFIS 
o 3 : ^RBF 
t ۱ 0 A OHolt-Winters 
E. | : OSARIMA 
Q0 ی یی‎ ۱ P ^ 
0 3 10 13 20 25 0 5 10 15 20 25 
25 


مقدار 


20 


"y 


^ 








TA 
Y 15 1 15 
3 10 
gt 
1 1 0۳15 
un a x OHolt-Winters 
2 - ۳ OSARIMA 
0 $ JA 0 ۰ 
0 5 10 15 20 25 0 5 10 15 20 25 
FDSI مقدار مشاهداتى شاخص‎ FDSI , مقدار مشاهداتی شاخص‎ 
25 

A 20 

1 15 

1 

33 10 

13 DANFIS 

ARBF 

z 5 po 

z Holt-Winters 

14 5ه 8 = 

0 € : 
0 5 10 15 20 25 


مقدار مشاهداتی شاخص FDSI‏ 


مقايسة عملکرد مدلهاى سرى زمانى SARIMA‏ و Holt-Winters‏ با روش‌های هوس مصنوعى در پیش‌بینی AAY e‏ 


مقايسة ميانكين polio‏ مشاهداتى و پیش‌بینی‌شده بر اساس آزمون 2 (جدول ۵) نشان مىدهد كه فرض صفر مبنى 
بر برابربودن ميانكين سرى polio‏ مشاهداتى و ييش بينىشدة شاخص FDSD‏ در هيج كدام از ایستگاه‌های مورد مطالعه 
بر اساس روش‌های ANFIS‏ و RBF‏ در سطح خطای Y‏ درصد و بر اساس مدلهاى سرى زمانی SARIMA‏ و Holt-‏ 
۹ در سطح خطای ۵ درصد رد نمی‌شود. بنابراين» مىتوان نتيجه كرفت كه don‏ روش‌های استفادهشده (بهويزه 
روش‌های ANFIS‏ و (RBF‏ به‌منظور پیش‌بینی شاخص FDSD‏ ميانكين سرى زمانی مشاهداتى را حفظ كردهاند (جدول 


۵ اين ویژگی می‌تواند تاکیدی بر کارایی روش‌های هوش مصنوعی در پیش‌بینی فراوانی روزهای همراه با طوفان 


گردوغبار باشد. 
جدول ۵. آزمون مقايسة میانکین سری زمانی مشاهداتی و پیش‌بینی‌شده شاخص FDSD‏ 
EUM ۲‏ آمارخ 2 
ایستگاه ‏ , ۳ Holt- Holt-‏ 
مشاهدا 5 ANFIS RBF SARIMA ANFIS 1 : SARIMA‏ 
Winters Winters‏ 
aha‏ د ۳۹۷ PA. PEE A‏ ۲ ۴ ۶ ۵ ۷ 
خاش ۴ ا وكيد “الاب EE WEE . mc‏ 
ایرانشهر YY ۰/۳۱۵ ۰۴١۱ EA «(NAA YYY ۰۴۵۸۵۵ EA UM‏ 
زاهدان ۵۰۲+ ARAN A YYY A ١ ۰۵١ UNE ٠ «(NAA‏ 
زابل «NY YAY «YYY ۰۶۹۸ +۵۲ «AM‏ ۰*۵۱ ۰/۶۱۵ ۰۷۰۵ 
نتيجه کیری 


هدف از اين پژوهش مقايسة عملکرد مدل‌های سری زمانی SARIMA‏ و Holt-Winters‏ با روش‌های هوش مصنوعی 
شامل شبكة عصبی مبتنی بر توابع ib‏ شعاعی (RBF)‏ و سیستم استنباط عصبی- فازی تطبیقی (ANFIS)‏ به‌منظور 
پیش‌بینی فراوانی روزهای همراه با طوفان‌های گردوغبار (FDSD)‏ در مقیاس فصلی در ينج ایستگاه سینوپتیک استان 
ks n‏ بكو زع و اه usu Ve EE‏ ایس فا al‏ کش یه اا EA‏ مود 
مطالعه روش ANFIS‏ بهترین عملکرد را نسبت به سایر روش‌ها داشت. همچنین. با افزايش شاخص FDSD‏ در 
ایستگاه‌های مورد مطالعه» دقت پیش‌بینی همه روش‌ها افزايش یافت. با اين توضیح که همه روش‌ها در ایستگاه زابل» که 
با ۷۱۱ روز در jb‏ زمانی ۲۵ساله (۲۰۱۴-۱۹۹۰) بیشترین فراوانی را به همراه داشت» بهترین عملکرد را از خود نشان 
دادند. غير از ایستگاه‌های خاش و سراوان» كه کمترین فراوانی روزهای همراه با طوفان گردوغبار را به‌ترتیب با مقادير 
۰ و ۱۰۶ روز در بازةٌ زمانی ۲۵ساله به خود اختصاص دادند. در دیگر ایستگاه‌های مورد مطالعه. روش‌های ۳15( 
Holt-Winters RBF‏ و SARIMA‏ به‌ترتیب رتبۂ اول تا چهارم را در پیش‌بینی شاخص FDSD‏ به‌دست آوردند. با کاهش 
فراوانی روزهای همراه با طوفان گردوغبار در ایستگاه‌های سراوان و خاش مدل سری زمانی Holt-Winters‏ به‌ترتیب 
عملکرد تقريباً یکسان و بالاتری نسبت به روش RBF‏ از خود نشان داد که حاکی از قابلیت بالای اين مدل به‌منظور 
پیش‌بینی مقادیر كم شاخص FDSD‏ بود. همچنین» نتایج بيانكر اين موضوع بود كه مدل سری زمانی SARIMA‏ نسبت 
به ساير روش‌های بيش بينى قابلیت بالایی در پیش‌بینی شاخص FDSD‏ در هیچ یک از ایستگاه‌های مورد مطالعه نداشت. 
همچنین, با وجود وقوع کمتر فراوانی روزهای همراه با طوفان گردوغبار در ایستگاه ایرانشهر نسبت به ایستگاه زاهدان» 
همه روش‌های پیش‌بینی شاخص FDSD‏ بر اساس معیارهای ارزیابی» عملکرد بهتر و دقیق‌تری نسبت به ایستگاه زاهدان 
داشتند که علت أن را می‌توان در وجود سری کامل و بدون انقطاع شاخص FDSD‏ در ایستگاه ایرانشهر جست‌وجو کرد؛ 


به طورى كه در ايستكاه ایرانشهر فقط در دو مقطع زمانى تابستان و jek‏ سال‌های ۲۰۰۵ و ۲۰۱۱ دو فصل متوالى عارى 
از روزهاى همراه با طوفان گردوغبار رخ داد. از طرفی» در پیش‌بینی شاخص FDSD‏ در استان سيستان و بلوجستان» مدل 
بهينة پیش‌بینی كننده مدل پیچیده‌ای بوده است. به صورتى كه برای همة ايستكادهاى مورد مطالعه. مدلى که از سه يا 
کار كام eb‏ هر ن ےم استفاده کرد هعتوان رین مدل فش كدو os‏ في از این رف دراك يدحا مائدة از 
طوفان‌های قبلی می‌تواند دلیل مهمی بر تأثير طوفان‌های چند فصل گذشته بر شکل‌گیری طوفان‌های گردوغبار در 
فصل‌های آینده باشد. يعنى طوفان شکل گرفته‌شده در حال حاضر ممکن است به خاطر ذرات رسوب کرده و ته‌نشین‌شده 
از سه يا چهار فصل گذشته باشد که اکنون به خیزش درآمده است. همین دلیل علمی که اکنون بیان شد در دیگر 
تاش ان رهق كاقة میا E hacc Sulis dafs‏ ان 5555 Bebidas.‏ 
سنجش صحیح عملکرد مدل‌ها در پیش‌بینی شده هم علت شکل‌گیری طوفان گردوغبار در مناطق مستعد اين پدیده را 
آشکار کرد؛ زيرا ذرات به‌جامانده از طوفان‌های گردوغبار قبلی خود عاملی بر شکل‌گیری طوفان‌های گردوغبار فعلی و آتی 


.+ 
هسنند. 
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